Magnetotransport and spin dynamics in an electron gas formed at oxide interfaces 
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We investigate the spin-dependent transport properties of a two-dimensional electron gas formed 
at oxides' interface in the presence of a magnetic field. We consider several scenarios for the oxides' 
properties, including oxides with co-linear or spiral magnetic and ferroelectric order. For spiral 
multiferroic oxides, the magnetoelectric coupling and the topology of the localized magnetic moments 
introduce additional, electric field controlled spin-orbit coupling that affects the magneto-oscillation 
of the current. An interplay of this spin-orbit coupling, the exchange field, and of the applied 
magnetic field results in a quantum, gate-controlled spin and charge Hall conductance. 

PACS numbers: 75.70.Cn, 73.40.-c, 72.25.-b, 68.47.Gh, 73.20.-r, 85.75.-d 



I. INTRODUCTION 



The transport properties of a semiconductor-based 
two-dimensional electron gas (2DEG) by external elec- 
tric and magnetic fields is at the heart of mesoscopic 
and spintronic research with a wide range of applica- 
tions. A new impetus has been the recent discovery 
of 2DEG formed at the interface of insulating oxides^. 
This 2DEG can be laterally confined and patterned to 
achieve new functionalities such as oxide-based field effect 
transistors^. In view of the remarkable phenomena ob- 
served in the conventional 2DEG under a magnetic field, 
including the Shubnikov-de Hass (SdH) effect* 5 -^, and the 
spin/charge Hall effect^ in the presence of spin-orbit in- 
teractions (SOI)™, it is timely to consider the properties 
of oxide-based 2DEG in a magnetic field. Generally, the 
relatively large effective mass m* together with a high 
carrier concentration N e at the oxide interfaces (for in- 
stance m*/m e 1=3 3.2 with m e being the free electron 
mass and N e « 5 - 9 x I0 13 cto~ 2 at the LaA10 3 /SrTi0 3 
interface^) imply a strong magnetic field for the SdH 
oscillation and for the quantized Hall conductance to be 
observable. In addition, oxide interfaces have a multi- 
tude of inherent properties such as multiferroicity and 
strong electronic correlation o 1 1 ' 12 . It is our aim here to 
inspect how such properties affect the magnetotransport 
in 2DEG. One of the exciting results is that, utilizing 
2DEG formed at multiferroic oxides (e.g., RMnOaii^) 
the magnetotransport can be modulated by both a mag- 
netic field and also with a small transverse electric field 
(~ XkV/cm). The topology of the local helical magnetic 
moments in multiferroics results in an induced effective 
SOI 3 - that is electrically-^ and/or magnetically 15 tunable. 
Consequently, we find a large magnetic and electric field- 
dependence of the quantum oscillations in the longitudi- 
nal conductance, as well as a spin/charge Hall conduc- 
tance. Furthermore, the cases of collinear order and/or 
weak on-site correlations are also considered. 



II. THEORETICAL FORMALISM 

We start from the general case by considering a 2DEG 
formed at the intcrsurface of a spiral multiferroic oxide 
(e.g., RMnOa) with a polar oxide (cf. FigflJ. In a homo- 
geneous, static magnetic field B = (0, 0, B z ) the single- 
particle 2DEG Hamiltonian reads 



H = hk + h ex , with 
2m* 



(1) 



Un r ■ cr + .g/i^B • cr. 



Here g is the gyromagnetic constant, [Ib is the Bohr mag- 
neton, and cr are the Pauli matrices. Un r is the exchange 
field determined by Coulomb repulsion and Hund's rule 
coupling at the oxides ionic sites, and can be described 
by an effective unit- vector field n r and a strength U . A 
useful parameterizations is 

n r = (sin# r , 0, cos# r ) 
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FIG. 1. (a) Schematic view of 2DEG at the interface of a 
multiferroic oxide with a polar oxide and in the presence of 
a magnetic field B. The spiral (x — z) plane is perpendicular 
to 2DEG. (b) The 2DEG energy levels as functions of the 
scaled spin-orbit interaction a with A/Jkd = 20. The spin- 
down (gray) and the spin-up (black) subbands are separated 
by the exchange energy. 
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with 



where 



i r = q m • r 



and q m = (q, 0,0) is the spin wave vector of the 
spira l 13 ' 14 . 

A B = (-B z y, 0,0) 

in the Landau gauge. The case of a collinear order and/or 
weak electronic correlations follow as special cases, as 
discussed below. Upon a unitary gauge transformation 

U g = exp(—i9 r a y /2) 

the local quantization axis becomes aligned with n r on 
the expense of introducing the topological vector poten- 
tial 



-iWlVUg = (-hqa y /2,0,0). 



The transformed kinetic energy 3 - reads 



h k 



1 



2m* 



eAf 



A 9 ] 2 



(hereafter, transformed quantities are marked by a tilde). 
A g depends only on the geometry of the local magnetiza- 
tion at the oxide and acts as a helicity q and momentum- 
dependent effective SOI ~ qP x a y . The transformed term 
h ex turns however diagonal, h ex = Acr z (after omitting a 
small periodic modulation on the spins 1 ^). The value of 
A depends on the system under study: 

i) For a strong local correlation compared to the external 
Zeeman field we find A « U. This is the case of spiral 
multiferriocs. 

ii) The case of a collinear n(r) (no ferroelectric order) 



corresponds to q —> 0, and hence A„ 



0. 



in) For a weak local correlation field U — >• and/or un- 
der a strong magnetic field we find A qhbB z , i.e. the 
conventional 2DEG in a magnetic field. 
The cases i)- Hi) correspond to different combinations 
of the oxides at which interface the 2DEG is formed. We 
will focus here on i) and ii). Predictions for the case 
Hi) are readily deduced, the physics however in this case 
resembles the well-known situation in a semi-conductor- 
based 2DEG (without SOI). 
For the transformed system we infer 

[Px,h k + h ex ] = 0. 

Introducing the bosonic operator 

a = [P x — eB z y — iP y ]/ \j2rahw (2) 

and scaling all energies in unit of the cyclotron energy 

huj = heB z /m* , 

we obtain the dimensionless harmonic Hamiltonian 

H = (a f a + 1/2) - aa y (a i + a) + a 2 + Ka z (3) 



A = A/fiw, a 



-, l B = \Jh/eB z . 



a is the lattice constant, and Ib is the magnetic length. 

In the limit of a collinear spin order, the eigenstates of 
the Hamiltonian are the Landau levels, 

\n,s) =^-L\0}®\s) 



with the energy 

e«s = (n+ 1/2) + sA + a 2 , n = 0, 1,2, • • • , s = ±. 

However, in most cases, the state \n, s) is coupled to 
\n ± 1, a) via SOI. The case of a weak exchange (local 
correlation) field follows for A = (cf. [l6j]). Expanding 
the eigenstate 

na 

in the Landau space, the secular equation 
leads to the following matrix form, 



(4) 







e n -i }S tasy/n 

iasy/n e n , s ias^Jn + 1 
ias\/n+ 1 6n+i,s 



The Hamiltonian in the Landau space reduces then to 
two independent, infinitely dimensional tridiagonal ma- 
trices with reference to two groups 17 \n, (— l) n ) and 
| n, — (— l) n ). The energy spectrum can be obtained nu- 
merically by truncating the matrix dimensions while in- 
cluding a sufficient number of Landau levels (up to 1000 
energy levels were taken during the calculations of the 
transport coefficients). 

FigHb shows the energy levels as a function of the SOI a 
with A = 20. Two spin subbands (parallel or antiparallel 
to the exchange field) are shown. Hence, the electronic 
structure varies strongly with a which is tunable by a 
small transverse electric field due to the magnetoelectric 
coupling 14 . The high density of states around a 5 fol- 
lows from a perturbation considerations: 
Let SOI be smaller than the exchange splitting, a« A. 
Up to a second order in a the two energy branches read 



E ns (a) 



At critical value 



(n + l)a 2 na 2 
s 1 — = — h s- 



2A 



[(4A 2 - 1)/4A] 



2A 



(5) 



1/2 



we find infinitely degenerate states, E ns= -i(a c ) = 0. For 
large a and/or n, the degeneracy is lifted and distributed 
around a c , e.g. a c = 4.47 with a given A = 20. 
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FIG. 2. The persistent spin current in units of J s o = 

vs. SOI strength a with A/hw = 43. Solid line: 
_Bf < A, only the spin-down subband intersects the Fermi 
level. Dotted line: Ef > A, both spin subbands are occupied 
when \a\ < 4. 



III. PERSISTENT SPIN CURRENT 

As pointed out in Ref.[3] without a magnetic field, a 
nonzero spin current (J^) is generated when only the 
spin-down subband is below the Fermi level Ep. With a 
magnetic field B z the spin current is 

(Jy) = iA/X>|#M/(J2*), (6) 




FIG. 3. The longitudinal conductivity a xx as functions of (a) 
magnetic field ujt and (b) SOI strength a at a fixed Fermi 
energy Eft /h = 60. The exchange energy is Ar/h = 43. 
Insets show the contributions from the different subband, 
is re-scaled with respected to the total conductivity a xx at (a) 
ujt = 0.5 and (b) a — 0. 



where v is the filling factor and / is the Fermi distribu- 
tion function. The spin current operator J| for the y 
component of the spin is defined as 



J v s = ^(.VyVx + v x a y ), 



[x,H] 

ih 



T. (O 

2to lv 



2acr„, 



(7) 
(8) 



Note, d y = dy because [a y ,U g ] = 0. An overall in- 
dependence of the amplitude of spin current is related 
to the electron density through the filling factor (Fig^J) ■ 
In the case of Ep/hu> = 60, two subbands are both occu- 
pied only when |a| < 4. It is then numerically confirmed 
that the persistent spin current vanishes when two spin 
subbands are intersected by the Fermi level^. No charge 
current is generated, i.e. ( J cx ) = (—ev x ) — in the ab- 
sence of an in-plane electric field. The spin polarization 
vanishes as well (a y ) = under a normal magnetic field; 
only a pure, electrically tunable persistent spin current 
exists. 



IV. LONGITUDINAL CONDUCTIVITY 

The topological SOI is exactly analogous to the 
semiconductor-based 2DEG with the Rashba and Dres- 
selhaus spin-orbit couplings being equal in strength. 



Therefore, the spin beats in the magneto-oscillation at 
the helimagnetic interface are completely suppressed 18 . 
However, the oscillatory magnetoresistance contains at 
least two components: SdH oscillations of the spin-down 
and spin-up subbands, respectively^. A 3> 1 means 
a large energy separation between the two subbands, 
~ 2(<7 Z )A, which results in negligible off-diagonal terms 
of the Green function^. The longitudinal conductivity 
is then given by a xx 



u XX 



Oxx with 



r(±) 



muj \ - 



{ll>±\JcxW±){ll>'±\Jcx\i>±) 



An 2 ^ {E^ ± - fi± + ih/r±)(E^ - fi± + ih/r±) 



where {^±\J C x\' l P'±) is the matrix elements of the charge 
current density operator. /i± is the energy distance be- 
tween the Fermi energy and the bottom of the spin- 
uptown) subband. t± is the total relaxation time in- 
cluding the intrasubband and intersubband scattering 
at a short-range random potential. Experiments on 
semiconductor-based 2DEG 21 suggest a Lorentzian shape 
as an excellent fit for the disorder-induced broadening 
of Laudau levels; we assume this to hold in our case, 
i.e. t± = t. Fig|3K demonstrates the dependence of 
<j xx on the magnetic field for the following parameters: 
Epr/h = 60, Ar/h = 43, and a = 0. The spectrum of 
a xx consists of two harmonics: The oscillating part of the 
partial conductivities a x ^ is described well analytically 
by cos(27T^- + 7r). The magnitudes are proportional to 
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the carrier concentrations in the spin-up and spin-down 
subbands, respectively. On the other hand, assuming a 
fixed Fermi energy the carrier concentration varies with 
the electrically tunable spin-orbit interaction parameter 
a, the quantum oscillation of the longitudinal conduc- 
tance a xx is caused then by a crossing of the chemical 
potential and the energy levels, see also in FigfJb. The 
harmonics stem from a change of the energy distance /x± , 
which induces an equivalent frequency of a X x ■ As a > 4, 
only the low subband is occupied however, ait^ — > 0. 

When accounting explicitly for intersubband transi- 
tions, the total relaxation time t t reads 2 - 

r- = r[l + (1 - + ^F+] (9) 

T T 

r+ = r[l + (1 - ^)F+ + - f F^ (10) 

T T 

with the times r and t' for intra- and inter-subband scat- 
tering at a short-range random potential, respectively. 
F± is oscillatory as inferred from 

F± = 2 cos(2tt^± + tt) exp( — — ) (11) 

nw LOT 

to the first order of the parameter exp(— — ). FigfJ] 
presents the dependencies of the a xx on the magnetic 
field under the condition that only the lower subband 
is below the Fermi level: Eft/H = and At/S = 43. 
A high-frequency oscillation component emerges (termed 
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FIG. 5. (Top): The charge Hall conductivity, a xy vs. SOI 
parameter a with no in-plane component of the magnetic field 
(By/Bz — 0). Inset: a close-up of the quantization of <j xy . 
(Bottom): The spin Hall conductivity, a s under B y /B z — 5. 
The exchange energy is A/hw = 43. 
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FIG. 4. The longitudinal conductivity a xx vs. magnetic field 
lot for various intersubband scattering intensities: (a) t/t' — 
0.0, (b) t/t' = 0.1 at a fixed Fermi energy E f t/H = 0. The 
exchange energy is Ar/ft = 43 and the SOI a = 0. 



magneto-intersubband scattering (MIS) effect^), which 
behaves as cos(4-7r ) and is distinguishable from the 
background harmonic cos(27r^j + tt). Upon increasing 
the temperature, MIS dominates the oscillatory magne- 
toresistance due to its weak-temperature damping com- 
pared to SdH oscillations^. 



V. HALL EFFECT 

Following Rashba 2 ^, the charge (spin) Hall conductiv- 
ity is given by Kubo-Greenwood formula, 



x y, s 



Til 2 ^ 
B 



E, 



</>' 



(12) 



The prime over the sum indicates a summation only over 
the state \ip') below the Fermi level (E^r < Ep) and 
above the Fermi level (E^ > Ep). The numerical 
evaluation shows a quantization of the charge Hall con- 
ductance, a xy — ve 2 /2-kH (see in FigfSh.), for the charge 
current carried by each state is not changed by SOI 7 . 
A similar electric-field effect on the charge Hall conduc- 
tivity has been reported in Graphene p-n junctions^, in 
which the carrier type and density can be controlled by 
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external gates as well. Different from a semiconductor 
2DEG, where a resonant spin Hall conductance is predi- 
cated for equal Rashba and Dresselhaus SOI 2 -, the spin 
Hall conductivity is zero in our system in the presence of 
only a perpendicular magnetic field due to the absence 
of a crossing of any adjacent energy-levels. 

When the magnetic field is tilted in the y — z plane 
with respected to the z axis, the in-plane y compo- 
nent By induces a nonzero spin polarization (<j y ) = 
^ / v^2^{%p\(7y\il}) , which renders the oxide magnetic or- 
der not exactly coplanar. A nonzero Berry's curvature 
is expected in the momentum spaced, and so is the spin 
Hall conductivity a s . The a-dependence of o~ s is dif- 
ferent from that of o~ xy (see Fig(S}}), which implies that 
the spin Hall conductance is not totally dominated by 
the spin-polarized charge current. A step- like structure 
is found in the spin-Hall conductivity for Ep < A. It 
should be noted that the filling factor v, the persistent 
spin current (Jf), and the charge Hall conductivity a xy 
are not sensitive to the in-plane magnetic field B y . How- 
ever, the competition between the Zeeman energy and 
the SOI results in sharp changes of the spin polarization 
(<T y ), which gives rise to resonant peaks around the de- 
generate point a c (for Ep = 0) and to a quick oscillations 
(for Ep — GOhaj) of the spin Hall conductivity a s . 



VI. SUMMARY 

In summary, we studied the spin current, the magneto- 
oscillations, and the spin/charge Hall effect in 2DEG at 
oxide interfaces exposed to magnetic fields. The trans- 
port behavior is determined by an interplay between the 
exchange interaction, the magnetic field, and the effective 
SOI induced by the spiral geometry of the local magnetic 
order. The electrically tunable topological SOI allows 
a control of the system's magnetotransport by a small 
transverse electric field. Our predictions are accessible 
experimentally. In fact, very recently, the SdH effect 
has been experimentally observed for SrTiOa/LaAlOa 27 , 
where an intense magnetic field up to 31. 5T was applied. 
A collinear spin phase is thus achieved, and the spin helic- 
ity q — > 0, which corresponds to a special case described 
by the Hamiltonian Eq(3]with a — and A = gfiBB z . 
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